The yielding dynamics of a colloidal gel 
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Attractive colloidal gels display a solid-to-fluid transition as shear stresses above the yield stress 
are applied. This shear-induced transition is involved in virtually any application of colloidal gels. 
It is also crucial for controlling material properties. Still, in spite of its ubiquity, the yielding 
transition is far from understood, mainly because rheological measurements are spatially averaged 
over the whole sample. Here, the instrumentation of creep and oscillatory shear experiments with 
high-frequency ultrasound opens new routes to observing the local dynamics of opaque attractive 
colloidal gels. The transition proceeds from the cell walls and heterogeneously fluidizes the whole 
sample with a characteristic time whose variations with applied stress suggest the existence of an 
energy barrier linked to the gelation process. The present results provide new test grounds for 
computer simulations and theoretical calculations in the attempt to better understand the yielding 
transition. The versatility of the technique should also allow extensive mesoscopic studies of rupture 
mechanisms in soft solids ranging from crystals to glassy materials. 



I. INTRODUCTION 

Colloidal particles with strong enough short-range at- 
traction are known to aggregate into fractal clusters that 
may form a space-spanning network leading to solid-like 
properties even at very low volume fractions [l|, Q ■ Such 
gelation obviously has tremendous consequences in terms 
of viscosity increase and rheological properties, which are 
capital for applications ranging from paints and coatings 
to food products and cosmetics. In the gel state, the 
system is kinctically arrested and displays non-ergodic 
features such as ultraslow relaxations, aging, and dy- 
namical heterogeneities 0, 0, H| ■ These glassy- like fea- 
tures have raised the issue of unifying the behaviour of 
colloidal gels and that of more concentrated amorphous 
jammed systems, e.g. colloidal glasses or concentrated 
emulsions [f|. Therefore, most recent works on systems 
of short-range attractive particles have focused on the 
physics ofgelation and on their structural properties at 
rest 00, HE! Ell- 

Of even greater practical importance is the behaviour 
of such materials under external stress. It is well known 
that applying stress to a soft glassy material generally 
triggers a solid-to-fluid transition, in the same way as 
increasing the temperature T, lowering the interparti- 
cle attractive energy U, or lowering the particle vol- 
ume fraction <f> [(I [12]. In particular, the yielding phe- 
nomenon, i.e. the fact that the material flows like a liq- 
uid above a yield stress a y and remains solid below a y , 
has been extensively documented in colloidal gels and 
glasses through standard rheological measurements and 
theoretical or numerical models, with emphasis on the 
scaling of viscoelastic properties and of Oj, with <f> and 
U 0,1, IH El, EI Eil E3, El El, H E3- However, 
rheology provides information that are averaged over the 



whole sample size and it cannot capture local phenomena 
such as apparent wall slip [221, sh ear localization [23| , or 
fracture within the material [2J]. Recently, a few orig- 
inal experimental approaches have been reported using 



confocal microscopy [25|, [2g. 
locimetry in microchannels 
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and particle imaging ve- 
. Still, these works were 
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restricted to steady-state measurements, and, as a mat- 
ter of fact, very little is known on how the stress-induced 
solid-to-fluid transition exactly proceeds in both space 
and time. 

In the present article, we introduce high-frequency ul- 
trasound as a new tool to follow the local deformation 
and flow of carbon black suspensions, a totally opaque 
colloidal gel made of weakly attractive soot particles. 
Both creep tests and oscillatory shear experiments show 
for the first time that, above a critical yield stress a y , (i) 
yielding proceeds from the cell walls in a spatially and 
temporally heterogeneous way and (ii) the time needed 
for total fluidization of the sample decreases exponen- 
tially with the applied shear stress. These results point 
out the existence of an energy barrier that linearly de- 
creases with the external stress. Therefore, they pro- 
vide evidence for the relevance of activated barrier hop- 
ping and for processes reminiscent of heterogeneous nu- 
clcation, in strong support of recent models [2911301]. Such 
features should be common to all weakly attractive col- 
loids and models are expected to become quantitatively 
predictive, based on more spatially and temporally re- 
solved measurements as described here. 



II. MATERIALS AND METHODS 

A. Carbon black gel preparation 

Carbon black is widely used in industry, e.g. for coat- 
ings, printing, rubbers, tires, paints, or batteries [3~j| . 
These colloidal carbon particles result from partial com- 
bustion of fuel and are made of unbreakable aggregates of 
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permanently fused nanomctric primary particles. These 
aggregates have a typical diameter of 500 nm. When dis- 
persed in a light mineral oil (from Sigma, density 0.838, 
viscosity 20 mPa.s) at a weight concentration of a few 
percents, carbon black particles (Cabot Vulcan XC72R 
of density 1.8) form a space-spanning network of fractal 
dimension df = 2.2 due to weakly attractive interactions 
of typical strength U ~ 30/csT Q. Our samples are 
prepared by vigorous mixing of 6 % w/w carbon black 
particles into the mineral oil. Since the resulting gel is 
transparent to ultrasound, 1 % w/w hollow glass micro- 
spheres of mean diameter 6 /im (Sphericel, Potters) are 
added in order to provide enough acoustic contrast (see 
below). The sample is further sonicated for 1 hour. 



B. Standard rheology 

Rheological measurements are performed in a Plexi- 
glas Couette cell (rotating inner cylinder radius 24 mm, 
gap width 1 mm, and height 30 mm) by a stress-imposed 
rhcomctcr (Bohlin C-VOR 150, Malvern Instruments). 
To ensure an initial reproducible gel state, all measure- 
ments are carried out by first preshearing the suspension 
at 1000 s _1 then at —1000 s^ 1 for 20 s each, in order to 
break up any large aggregate. The sample is then left at 
rest for 100 s during which the gel structure reforms, a 
process that is controlled by monitoring the viscoelastic 
moduli G' and G" at / = 1 Hz and a = 0.2 Pa. A rest 
time of 100 s is sufficient for G' and G" to reach a steady 
state with negligible aging (Fig. [2ji). All measurements 
were performed at a temperature of 25°C. 



C. Ultrasonic echography under shear 

In our experiments, classical rheology in Couette ge- 
ometry is combined to ultrasonic echography in order to 
access local properties of the gel such as its velocity or de- 
formation fields. In creep experiments, we use ultrasonic 
speckle vclocimetry (USV) as described in Ref. [12] to 
record temporally-resolved velocity profiles accross the 
gap at about 15 mm from the cell bottom (Fig. QJ,). 
USV is based on the interaction between ultrasound and 
micronsized scatterers embedded in the gel (here, the 
glass microspheres that act as contrast agents). A high- 
frequency (36 MHz) immersion transducer first emits a 
short ultrasonic pulse and collects the pressure signal 
backscattered from the gel. The backscattered signal is 
an ultrasonic speckle that directly reflects the spatial dis- 
tribution of the scatterers within the gel along the acous- 
tic beam (Fig. [I})). As the scatterers move with the gel, 
cross-correlating two successive speckle signals over small 
time windows gives access to the local velocity of the gel 
with a spatial resolution of about 40 /im. This original 
combination of rheology and ultrasound turns out to be 
very efficient to follow the onset of flow both spatially 
and temporally. 
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FIG. 1: Methods: ultrasonic echography under shear, (a) 
Schematic of the rheology experiment combined with USV. 
Temperature is controlled to 25±0.1°C by a water circulation 
around the Couette cell, (b) Typical normalized pressure sig- 
nal backscattered by glass microspheres embedded in a car- 
bon black gel as a function of time after the incident pulse is 
emitted at t = s. The times-of-flight corresponding to the 
positions of the stator and the rotor are indicated by a red 
and blue dashed line respectively. The corresponding gap is 
1 mm. (c and d) Spatiotemporal diagrams of the normalized 
pressure signal during an oscillation experiment at / = 1 Hz 
and a = 15 Pa. The pressure is coded in linear grey levels and 
plotted against the time at which the incident pulse is sent 
(x-axis) and the radial distance r to the rotor (y-axis). In (c), 
the pulse repetition frequency / pr f is 400 Hz, which is large 
enough to resolve the deformation due to the oscillations of 
the rotor, while in (d), / pr f = / = 1 Hz. (e and f). Same as 
in (c and d) but for a higher shear stress a = 30 Pa. 
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For oscillation experiments, we take advantage of ul- 
trasonic echography in a new way by focusing on the 
spatiotcmporal diagrams of the successive speckle signals 
(Fig. [Iji-f). Indeed, the spatial variation of the speckle 
intensity provides crucial information on the local defor- 
mation of the sample. When pulses are sent with a suf- 
ficiently high repetition frequency / pr f compared to the 
oscillation frequency /, tracking the local oscillatory mo- 
tion of the gel becomes possible (Fig.QJ and e). Here, we 
rather set / pr f = / so that the material is probed by ultra- 
sound only once per oscillation period. Data analysis is 
then straightforward: if two successive speckle signals are 
fully correlated, then the gel is solid-like since it comes 
back exactly to the same position every period (Fig. [TJi)- 
On the other hand, we interpret any significant degree of 
decorrelation in the speckle signal from one period to the 
other as the signature of fluid-like behaviour due to mo- 
tions of the scatterers induced by local rearrangements 
inside the gel (Fig. [If) . Note that such decorrelation in 
the fluidized state does not arise from flow irreversibil- 
ity (since the Reynolds number is Re ~ 10 -3 <C 1) but 
rather from sedimentation of carbon black clusters due 
to the density mismatch between the mineral oil and 
the carbon black particles. We checked that sedimen- 
tation of clusters as small as 10 /im is sufficient to pro- 
duce significant decorrelation after one oscillation period. 
The spatiotemporal tracking of the transition from corre- 
lated speckle signals (i.e. solid-like state) to decorrelated 
speckle signals (i.e. fluidized state) is made quantitative 
by using a threshold for the correlation coefficient of suc- 
cessive signals over small time windows, thus identifying 
the boundary between solid and fluid regions as a func- 
tion of time up to a spatial resolution of 40 /im (Fig. [Hi). 

The fraction of fluidized sample <f> at a given time t 
is readily extracted from the red contour separating fluid 
and solid regions in Fig.rjji. Jumps in <j)(t) are interpreted 
as fast fluidization events and identified using a threshold 
on the second derivative of <j>(t). Statistical analysis of 
the time intervals St between two successive jumps then 
yields the probability distribution function P and the 
complementary cumulative distribution function C of St 
(inset of Fig. [5]). C is defined as C(x) = f P(y)dy and 
experimentally estimated by sorting the measured values 
{Stj}j = i...N in ascending order and by plotting 1 — j/N 
as a function of Stj I33II. 



III. RESULTS 

A. Standard rheology 

Our experiments focus on a colloidal gel made of car- 
bon black particles dispersed at 6 % w/w in a mineral oil 
[HI . Its rheological properties are typical of a soft solid as 
shown in Fig.[2j the elastic modulus G' is about ten times 
larger than the loss modulus G" and both moduli are fre- 
quency independent (Fig. [5^). Moreover, as the ampli- 
tude a of the oscillatory shear stress is increased, the gel 
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FIG. 2: Standard rheology of a 6 % w/w carbon black suspen- 
sion in a mineral oil. Viscoelastic moduli G' (■) and G" (O) 

(a) in a frequency sweep in the linear regime at a = 0.5 Pa, 

(b) in a stress sweep experiment at / = 1 Hz and (c, d) at 
a — 0.2 Pa and / = 3 Hz as a function of time after preshear 
at high speed during 60 s. The red dashed line in (b) indicates 
the yield stress a y , which is time and frequency independent. 



yields for a = a y ~ 8 Pa (Fig.rjb), where a y is defined by 
G'(a y ) = G"{(7 y ). Above cr y , one has G" > G', i.e. the 
material behaves like a liquid. As seen in Fig. [2t where 
the gel was fluidized by strong prcshcaring for times t < 
and left at rest for t > 0, the elastic modulus recovers 
within a few seconds. This variation corresponds to the 
gelation process. The short range attraction between the 
carbon black particles leads to cluster aggregation. As 
the clusters percolate, they form a space-spanning net- 
work. The backbone of the network gives rise to the 
observed solid-like structure. For t > 100 s, we observe 
very little aging (Fig. [2U), which is ideal to run long ex- 
periments. 

Yield stress measurements and global rhcological prop- 
erties are, however, too simple and elude fundamental is- 
sues such as the mechanism of fluidization and the role of 
self-healing processes due to the attractive interparticle 
forces. To follow the local dynamics of yielding, rheome- 
try in Couette (concentric cylinders) geometry is coupled 
to ultrasonic speckle velocimetry (USV, see Methods sec- 
tion), which allows us to access deformation and velocity 
fields with high temporal resolution [32j|, under either 
creep tests or oscillatory shear. 
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B. Creep experiments 

Figure^ gathers the results of creep experiments for 
a wide range of imposed shear stresses. The critical shear 
stress <x y = 8.5 ± 0.5 Pa clearly separates two regimes. 
For a < a yi the gel remains solid-like and does not flow, 
whereas for a > a y , the sample eventually flows with a 
measurable shear rate 7. The shear rate response j(t) 
presents a characteristic shape that has already been 
observed in creep experiments on surfactant columnar 
phases [35| . In order to get better insight into the flu- 
idization process, we correlate the rheological measure- 
ments to the local velocity field v(r) measured with USV 
for the typical creep experiment at a = 10 Pa shown in 
Fig. [3Jd. In the early stage where j(t) slowly increases 
(t < 1700 s), solid-body rotation is observed, i.e. the 
sample is in an unsheared solid-like state except for thin 
lubrification layers along the cell walls that allow the gel 
to slide as a solid block (Fig. [3):). At the beginning of 
the sharp upward bend in -y(t), for 1700 s < t < 1800 s, 
the bulk sample starts to experience a non-zero shear 
rate. Fluidization then proceeds in a spatially hetero- 
geneous way. Indeed, in this regime, velocity profiles 
present a highly sheared region close to the rotor that 
coexists with a low-shear band next to the stator and 
whose size increases in time (Fig. [3jl). Such shear local- 
ization suggests that the gel progressively sticks to the 
cell walls and gets dragged along the rotor. The fluidized 
zone expands from the rotor until, for t > 1800 s, the gel 
flows homogeneously without any apparent wall slip, as 
seen from the linear velocity profiles of Fig. [3p. In this 
last stage, the increase of the bulk shear rate indicates 
that further microstructural changes still occur and/or 
that sedimentation sets in. 

Our local measurements refute the simple scenario in 
which yielding occurs abruptly and homogeneously at a y 
(as suggested by the standard stress sweep experiment of 
Fig. ^Bp) and rather support a yielding mechanism that 
is heterogeneous in both space and time, where fluidiza- 
tion starts at the walls and propagates through the whole 
sample. Such front dynamics is reminiscent of hetero- 
geneous nucleation as predicted by a phenomenological 
"fluidity" model (3|| . The present experiment also allows 
us to define unambiguously a characteristic time r for 
fluidization as the time when homogeneous shear is first 
observed, e.g. r = 1850 ± 100 s for a = 10 Pa. The flu- 
idization time roughly corresponds to the inflection point 
in the 7(4) curve (and typically to 7 = 1-10 s _1 ), a crite- 
rion that has been used previously [35| but without the 
confirmation from local measurements. 

Finally, Fig. [3f sums up the behavior of our colloidal 
gel in creep experiments. Below a y (region I), the gel 
remains solid for all times, while above a y , the gel first 
undergoes solid body rotation or shear localization for 
t < t (region II) and eventually flows homogeneously 
for t > t (region III). Remarkably, the fluidization time 
follows an exponential law: r ~ cxp(— u/cfq) where <To = 
4.0±0.2 Pa. In other words, with ctq = JcbT/v, where v is 
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FIG. 3: Creep experiments in a 6 % w/w carbon black gel. (a) 
Temporal evolution of the shear rate 7 for (from bottom to 
top) a = 5,7, 8, 9, 10, 11, 13, 15, 18, 20, 23, 25, 28, 30, 33, and 
40 Pa. (b-e) Combined rheological and USV measurements 
for a — 10 Pa. (b) 7(t) as a function of time. The red dashed 
line indicates the fluidization time r discussed in the text, (c, 
d, e) Instantaneous velocity profiles v(r) recorded with USV 
at the times shown in (b) using the same symbols, (f) The 
fluidization time r as a function of the applied shear stress 
a. The black line is an exponential fit to the experimental 
data. Region I: solid-like behaviour. Region II: plug-like flow 
or shear localization. Region III: homogeneous flow. 
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some characteristic volume, this behaviour is equivalent 
to an Arrhenius law r ~ exp(E(a)/(ksT)), where the 
energy barrier E(a) = —av decreases linearly with the 
applied stress. 



C. Oscillatory shear experiments 

In a second series of experiments, the gel is submit- 
ted to oscillatory shear stress of given amplitude a and 
frequency / = 1 Hz. The strain amplitude 70 (t) is moni- 
tored and again, we observe that when a is large enough 
70 eventually increases with time, suggesting that the 
sample becomes fluid over time (Fig. |4ji,b). To confirm 
this scenario, we track the local deformation of the gel 
with USV (see Methods section). For a < 5 Pa, the gel 
indeed remains solid: no decorrclation of the ultrasonic 
speckle signals, i.e. no rearrangement within the gel, was 
observed during waiting times as long as 10 5 s. How- 
ever, for stress amplitudes above 5 Pa, rearrangements 
and fluid-like behaviour are detected first close to the 
cell walls (see t ~ 100-300 s in Fig. Hi). The fiuidized 
zone then progressively expands through the bulk ma- 
terial (t ~ 300-550 s) until the whole sample gets fluid 
(t > 550 s). 

It is important to note that fluidization does not have 
such a dramatic signature on 70 (i) under oscillatory 
stress as on j(t) during creep tests. Moreover, as shown 
in Fig.[3h, the loss modulus is larger than the elastic mod- 
ulus throughout the whole experiment. In view of the 
USV results, this can only be explained by the presence 
of lubrication layers at the walls. A naive interpretation 
of standard viscoelastic data would thus wrongly con- 
clude that the sample is fluid as soon as stress is applied. 
Here simultaneous spatially-resolved measurements show 
that solid-fluid coexistence comes into play and reveal a 
highly heterogeneous fluidization pattern. 

In oscillation experiments, the fluidization time, de- 
fined as the time r for which 90 % of the gel is in the fluid- 
like state, follows an exponential law r ~ exp(— u/ctq) 
with cr = 3.9 ± 0.3 Pa for a > 8 Pa (Fig. gt). This 
is strikingly similar to the behaviour obtained in creep 
experiments (Fig. |3f). However, under an oscillatory 
stress of amplitude 5 Pa < a < 8 Pa, a strong depar- 
ture from exponential behaviour is observed, pointing to 
a divergence at a ~ 5.5 Pa. In this stress range, no flu- 
idization was observed during creep tests. Although this 
regime certainly deserves more attention, in the follow- 
ing, we shall focus on the exponential behaviour observed 
for stresses above 8 Pa. 

Shown in Fig. [5] is the fraction of fiuidized sample 
4>{t) for ten different oscillatory experiments at a fixed 
stress amplitude of a = 15 Pa. The large deviations 
from one run to the other reveal the random nature 
of the local processes at work during yielding. These 
noisy features are compatible with models that include 
activated dynamics [l6|, H3, H3- Overall, the fiuidized 
fraction (4>(t)) averaged over the various experiments 
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FIG. 4: Large amplitude oscillatory shear experiments at / = 
1 Hz in a 6 % w/w carbon black gel. (a) Temporal evolution 
of the deformation amplitude 70 for (from bottom to top) 
a = 5, 5.8, 6, 7, 8, 10, 13, 15, 18, 20, 23, 26 Pa. (b-d) 
Combined rheological and USV measurements for a = 15 Pa. 
(b) 70 and (c) G' (•) and G" (•) as a function of time, (d) 
Spatiotemporal diagram of the ultrasonic speckle signal coded 
in linear grey levels. The USV sampling frequency is equal 
to the oscillation frequency / = 1 Hz. Red dots indicate the 
boundary between solid-like and fluid-like regions. The red 
dashed line shows the time r at which 90 % of the gel across 
the gap is fiuidized. (e) The fluidization time r as a function 
of the stress amplitude a. The black line is an exponential fit 
to the experimental data (only data above 8 Pa are considered 
in the fit). Region I: solid-like behaviour. Region II: solid- 
fluid coexistence. Region III: fluid-like response. Between 
the green and blue dashed lines, r strongly deviates from the 
exponential behaviour. 



is well fitted by a sigmoidal relaxation (4>)(t) = 1 — 
cxp(— £ 2 /t 2 ), a behaviour typically found in nucleation 
and one-dimensional growth phenomena [38| . 
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FIG. 5: Temporal evolution of the fluidized sample fraction 
<j> for ten different oscillatory shear experiments at / = 1 Hz 
and a = 15 Pa. Black dots show the average behaviour (4>(t)) 
and the thick green line is {<j>(t)) — 1 — exp(— t 2 /f 2 ) with 
f = 320 s. The data are plotted against t/r where r = 
480 s is the fiuidization time defined as (</>(t)} = 0.9. Insets: 
probability distribution function P (top) and complementary 
cumulative distribution function C (bottom) of the waiting 
times 8t between two successive jumps in cj>(t). The solid 
lines are exponential fits P(5t) ~ C(5t) ~ exp(— t/r*) with 
t* = 0.137r = 65.8 s. 



IV. DISCUSSION 

The creep tests and the oscillation experiments above 
have shown the benefit of local measurements in estab- 
lishing a spatiotemporal scenario for yielding: fiuidiza- 
tion is initiated at the walls and propagates into the en- 
tire sample within a characteristic time r that follows 
an Arrhcnius law with an energy barrier E = —<jv. It 
is interesting to notice that a similar Arrhenius scal- 
ing has been reported in shear-induced aggregation of 
non-Brownian suspensions [39| . In both cases of shear- 
induced fiuidization and aggregation, the energy barrier 
is reduced proportionally to the shear stress. 

To the best of our knowledge, the present results pro- 
vide the first experimental evidence for the relevance of 
activated processes and barrier hopping in the context of 
yielding in colloidal gels. Although such ingredients have 
recently been incorporated into theoretical approaches 
based on mode-coupling theory [l6[ , local dynamical de- 
scriptions are still missing from microscopic models of 
the stress-induced solid-to-fiuid transition. Our results 
also suggest a comparison with fracture in heterogeneous 



materials. Indeed, the propagation of fracture fronts is 
known to involve activated processes attibuted to the nu- 
cleation of microcracks [lo, ST], H^j. Similarly, yielding 
proceeds by rapid local fiuidization events (Fig. |3Ji) sep- 
arated by some "waiting times" St. Previous measure- 
ments on fibrous materials have shown that time inter- 
vals between microcrack events arc power-law distributed 
[4l| . a feature typical of self-similar avalanche- like phe- 
nomena. Here, St rather follows an exponential distri- 
bution (insets of Fig. [5]), which points to random mi- 
correlated events. Moreover, our observation of a linear 
energy barrier E = —ov is in stark contrast with the 
Griffith scaling E ~ a~ 4 for fracture [42j]. Following 
Ref. [1^, the linear dependence E = —ov may be in- 
terpreted as the reduction of the energy barrier through 
hydrodynamic energy. In that case, the volume v reads 
V = a 2 S with a the particle radius and 8 a microscopic 
length characteristic of the interparticle potential. Here, 
assuming that T can be identified with the thermody- 
namic temperature, ksT/v ~ 4 Pa and a ~ 500 nm lead 
to S ~ 4 nm, which provides a reasonable estimate for the 
range of attractive forces between carbon black particles. 



Very recently, a novel rheological model to describe 
the local dynamics of ductile to brittle fracture was in- 
troduced and applied in the context of metallic glass for- 
mers [43| , but the connection with soft gels remains to be 
made. In particular, it can be noted that in some cases 
the fluidized fraction <p(t) decreases over short time peri- 
ods before increasing again at longer times (Fig. [5]). This 
behaviour, absent from fracture phenomena, is a direct 
illustration of the competition between shear-activated 
cluster breaking (which favors fiuidization and an in- 
crease of <p(t)) and self-healing processes (which tend to 
lower 4>(t)) due to the attractive interactions. 



In summary, the fiuidization scenario depicted in this 
article highlights the crucial importance of timc-rcsolved 
local measurements in characterizing the behaviour of 
soft materials under stress. Ultrasonic echography con- 
stitutes a versatile tool to follow yielding and fiuidization 
in both space and time even in optically opaque materi- 
als. Its mesoscopic resolution of a few 10 ^m bridges the 
gap between microscopic imaging techniques and global 
rheological measurements. Indeed, coarse-grained mod- 
els for soft glassy materials use local strains and stresses 
averaged over intermediate sizes [3(1 H3] , so that gath- 
ering information at a mesoscopic scale is a key issue for 
comparing experiments to such models. From the fun- 
damental point of view, our observations not only em- 
phasize the importance of activated processes and noise 
in yielding but also reveal a complex interplay between 
boundaries and bulk dynamics. Such results should pro- 
vide a new test ground for computer simulations and for 
theoretical calculations in the current attempt to under- 
stand yielding in a wide variety of soft glassy materials. 
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